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The combination of plasmonic nanoparticles and graphene enhances the responsivity and spectral 
selectivity of graphene-based photodetectors. However, the small area of the metal-graphene junc¬ 
tion, where the induced electron-hole pairs separate, limits the photoactive region to sub-micron 
length scales. Here, we couple graphene with a plasmonic grating and exploit the resulting surface 
plasmon polaritons to deliver the collected photons to the junction region of a metal-graphene- 
metal photodetector. This results into a 400% enhancement of responsivity and a 1000% increase 
in photoactive length, combined with tunable spectral selectivity. The interference between surface 
plasmon polaritons and the incident wave introduces new functionalities, such as light flux attrac¬ 
tion or repulsion from the contact edges, enabling the tailored design of the photodetector’s spectral 
response. This architecture can also be used for surface plasmon bio-sensing with direct-electric- 
readout, eliminating the need of complicated optics. 


Graphene-based photodetectors (PDs)[I|,i have been 
reported with ultra-fast operating speeds (up to 262GHz 
from the measured intrinsic response time of graphene 
carriers i) and broadband operation from the visible 
and infrared (sl-ll^ up to the THz[iMi. The sim¬ 
plest graphene-based photodetection scheme relies on the 
metal-graphene-metal (MGM) architecture!^ 0, H, fill . 
[2 qI - [^ . where the photoresponse is due to a combina¬ 
tion of photo-thermoelectric and photovoltaic effects [HI, 
0,i im. Eor both mechanisms, the presence of a 

junction is required to spatially separate excited electron- 
hole (e-h) pairs [HI. ItI. IsI. fill. I2 qI- 221. At the metal-graphene 


junction, a work-function difference causes charge trans¬ 
fer and a shift of the graphene Eermi level underneath the 
contact y, [bl, [7|, l2^. cormared to that of graphene away 
from the contact 0, ialil, resulting into a build-m 
of an internal electric field (photovoltaic mechanism) [5|, 
and into a difference of Seebeck coefficients 
(photo-thermoelectric mechanism) [Tl|, [2l|, [2^ . An alter¬ 
native way to create a junction is to use a set of gate 
electrodes to electrostatically dope graphene [si flU. 


Eor both photovoltaic and photo-thermoelectric mech¬ 
anisms, however, the spatial extend of the junction 
is^l00-200nm[7l Isl, UJj, which reduces the photoactive 
area to a fraction of the diffraction limited laser spot 
size in a typical scanning current microscopy experiment. 
Eurthermore, suspended undoped graphene only absorbs 
2.3% of light [13 which, while remarkably high for a one 
atom thick material, is low in absolute terms for practi¬ 
cal applications. This is further reduced by a factor of 
4/(1+ n)^ for graphene on a dielectric substrate of refrac¬ 
tive index n (see Methods). Additionally, in highly doped 
graphene the absorption decreases even further |28l l29|. 


One approach is to extend the junction region in order 
to capture more light. In a vertical (i.e. with doping 
profile perpendicular to the device’s surface) p-i-n semi¬ 
conductor PD this is achieved by ion-implantation with 
tailored dose and energy [soj. In the MGM configuration. 


however, the lateral nature (i.e. doping profile parallel 
to the device’s surface) of the junction does not allow 
a straightforward doping profile engineering and, thus 
far, to the best of our knowledge, no techniques have 
been reported to reliably do that. As an alternative, sev¬ 
eral graphene-based vertical architectures have been pro¬ 
posed, including all-graphene [Hlj . or graphene integrated 
with semiconductor layers, such as other two dimensional 
(2d) materials [sii l33| or SiQ [HHI- In the latter cases, 
however, graphene is not the absorbing material and the 
spectral response is thus far limited to above the band 
gap of the semiconductor layer. Eurthermore, while these 
approaches have demonstrated high responsivities (up 
to^ 5 X lO^A/W in Ref. [13 by employing M0S2 ns light 
capturing material) they do come with the cost of smaller 
operation speed (up to lOOkHz in Ref. [HI) as compared 
to the graphene-based PDs operating at speeds up to 50 
GBit/s in the optical link reported in Ref. [13 • 

Improving graphene absorption in the ultra-fast MGM 
configuration is thus critical. Various solutions have 
been proposed, such as the integration of graphene into 
an optical micro cavity [13 113 (> 20-fold enhancement) 
or onto a planar photonic crystal cavity [13 (8-fold en¬ 
hancement), to take advantage of the multiple passes of 
the trapped light through graphene, or its coplanar in¬ 
tegration with a Si integrated photonic waveguide [iol - li^ 
(>I0-fold enhancement). Another solution is the integra¬ 
tion of plasmonic nanostructures on graphene [ol, ll0|, |^ 
to exploit the strongly enhanced electromagnetic near¬ 
fields |^, 1^ associated with the localized surface plas¬ 
mon resonances (LSPR)!!!, 113. LSPRs originate from 
the resonant coherent oscillation of the metal’s con¬ 
duction electrons in response to the incident radia¬ 
tion. The resulting enhanced near-fields surrounding the 
nanostructures promote light absorption in the materials 
around them [^. We previously reported a x 20 enhance¬ 
ment in photoresponse [3 when radiation is focused close 
to the nanostructures. In this approach, however, light 
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absorbed around nanostructures far from the junction 
(where efficient charge separation occurs) does not fully 
contribute to the photoresponse [i^. 

An ideal alternative would be to enable light collection 
in one part of a device and then guide it into the junc¬ 
tion region at the contact edge. This can be achieved 
by exciting surface plasmon polaritons (SPPs) on the 
metal contacts. SPPs are surface-bound waves propagat¬ 
ing on a metal-dielectric interface and originate from the 
coupling of light with the metal’s free electrons ^ M- 
Their excitation can be achieved by means of an inte¬ 
grated diffraction grating [iii li^, and their delivery to 
the active region (junction at contact edge) will enhance 
the overall absorption. Thus, the contact now becomes a 
light collector. Such an approach was demonstrated with 
semiconductor-based near infra red (NIR) PDs ^-53| in 
order to reduce the semiconductor active area without 
compromising light absorption. A smaller active area re¬ 


sults into reduced carrier transit time 48, ^ and reduced 


capacitance jdsi 1^. thus increased operation speed. In 
particular, a circular (bull’s eye) grating 48, was used 
to deliver SPPs into a subwavelength circular aperture 
on top of vertical SiQ or Ge[4^ Schottky photodiodes, 
while a linear grating was used to deliver SPPs into a 
sub wavelength linear slit in a lateral metal-GaAs-metal 
photodiode configuration [52[ [^. Op eration speeds be¬ 
yond lOOGHz where estimated |48l. l49| . while responsivity 
enhancements (compared to a device without the grat¬ 
ing) were up to x4 for linear gratings and over 

xlO for circular gratings |48[ l49| . No compromise in op¬ 
eration ^eed was reported due to the presence of the 
grating 
detection 


This approach was also applied in mid-IR 
^ ^ , by delivering light into a quantum cascade 
detector, and in THz detectionjHlI, by delivering light 
into a GaAs/AlGaAs 2d-electron-gas bolometer. 

Here we apply the SPP grating coupler concept to a 
MGM PD and demonstrate a^400% increase in respon¬ 
sivity and a^l000% increase in photoactive length. Fur¬ 
thermore, we show that this offers a solution to another 
problem: in order to have a net response under uniform 
illumination of the whole MGM PD area, one must break 
the mirror symmetry between the two contacts In 
contrast to the metal-semiconductor-metal case, apply¬ 
ing a bias is not practical because it would result into 
a large dark current, due to the semimetallic nature of 
graphene 0. Using different metallizations for the two 
contacts is an option^, but increases the fabrication 
steps. In our approach this problem is addressed by us¬ 
ing different contact grating structures. One can utilize 
the interference between SPPs and incident waves and 
create novel asymmetric contact designs that produce 
complex spectral responses, such as switching the light 
flux between the two contacts edges, enabling new func¬ 
tionalities, such as label- and optics-free direct-electrical- 
readout plasmonic biosensing. 

To explore the design opportunities offered by a SPP 


grating coupler on the metal contacts, we first carry out 
numerical simulations using the finite-difference time- 
domain (FDTD) method 13 (sgI. The SPP wavevector 
on a metal-dielectric interface is|45|: 

= (^/c)\/(e^^i^e^)7(D^r^r^ (1) 


where e^(co’), ed{uj) are the dielectric functions for the 
metal and dielectric medium respectivel y, a nd the SPP 
existence condition is HGj - kspp is 

larger than any propagating wave in the dielectric, whose 
wavevector is k = {uj/c)y^. This momentum mismatch 
between SPPs and propagating waves implies that SPPs 
cannot decay into free propagating waves, but also that 
one cannot directly excite SPPs from free waves on a 
smooth metal surface [iHl- Oiie way to overcome this is 
by diffraction, whereby the continuity of the component 
of momentum parallel to the surface is broken and inci¬ 
dent light can scatter into SPPs [3 • This can be achieved 
by a nano-slit 57^ 5^. a diffracting element (H^. or a 
grating coupler [59, 61|. In the latter case, in particular, 
a periodic array of metal ridges and grooves delivers the 
additional momentum according to|45j: 


kspp = ^11 A tuK (2) 

where = {uj/c) sin0 is the parallel component of in¬ 
cident wavevector, 0 the incident angle, K = 27r/a the 
grating’s reciprocal lattice vector, a the grating pitch, 
and m the diffraction order. 

We select a grating of 50nm Au bars periodically 
placed at a pitch of 620nm on top of a 50nm Au con¬ 
tact film, as depicted in FigH^, at a 1:1 ratio of ridge 
and groove widths. These grating parameters are chosen 
since they were shown to be optimal for yielding a high 
percentage of incident light scattered into SPPs (^20% 
in Refs. |591. Ibljb A termination ’’step” of width d ex¬ 
tends beyond the last ridge. For simplicity, we assume 
the graphene/contact structure to be on top of a semi¬ 
infinite Si02 substrate. The dielectric functions of Au[^ 
and graphene [3 treated through a Drude-Lorentz 
model, as explained in Ref. Q and in Methods. Insert¬ 
ing Au’s dielectric function in Eqs J1I2[ with ed=l for air, 
yields Ao=645nm for the vacuum wavelength of the SPPs 
on the Au/air interface. It is also possible to have SPPs in 
the Au/Si02 interface. Given that ed=2.13 for Si09[63|- 
they are excited at Ao=930nm. 

In the first set of calculations the device is illuminated 
by a normally incident plane wave polarized perpendicu¬ 
lar to the grating (transverse magnetic-TM). This polar¬ 
ization is required because the SPPs are themselves TM 
waves [3 liSl? having both longitudinal (parallel to the 
propagation direction) and transverse (perpendicular to 
the surface) electric field components. The exposed part 
of the single layer graphene (SLG) in-between the con¬ 
tacts is kept fixed at lOOOnm, while the width of the 
Au contact, thus of the SLG unexposed (buried) part, is 
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FIG. 1: (a) Schematic of the simulated structure, (b) Relative 
(i.e. normalized to the case without the grating) absorption in 
the exposed SLG. (c) Gumulative absorption in the exposed 
SLG integrated away from the contact edge, for the two cases 
pointed by the arrows. 


varied depending on the number of ridges and the size d 
of the termination step. Periodic and absorbing bound¬ 
aries are assumed in the lateral and vertical directions 
respectively (see Methods). FigdjD plots the absorption 
in the exposed SLG normalized to the flux incident on the 
exposed SLG area, for d=0 and 410nm. Strong enhance¬ 
ment peaks, reaching up to 5% absorption, are found be¬ 
tween 630-700nm. There are also secondary small peaks 
at^900-950nm, assigned to the Au/Si02 interface SPPs. 
The green dashed line indicates the absorption within 
the lOOOnm wide SLG on top of the semi-infinite Si02 
substrate in the absence of the grating, as derived in the 
thin film limit [see Methods]: 


Aslg = 4^Islg/ insi02 + 1)^ (3) 


Si 02 refractive index in the above wavelength range j^. 
A three-fold wavelength-selective increase in absorption 
is observed, due to coupling with SPPs scattered from 
the grating. 

We now consider the parameter d. FigdlD indicates 
that for certain wavelengths (e.g.^630nm), there are op¬ 
posite extremes of absorption for the two different steps. 
Such an asymmetry can be instrumental in designing con¬ 
tact layouts that exhibit a net photovoltage even under 
uniform illumination. For the two cases pointed by the 
two arrows of Fig (Ud, we report in graph Fig (It the cumu¬ 
lative absorption in the exposed SLG as we move away 
from the grating edge. In the d=410nm case, absorption 
is strongly enhanced at the edge of the grating, and starts 
leveling off^300nm away from it. In contrast, for d=0 ab¬ 
sorption is minimal close to the grating and starts picking 
up only>300nm from the grating edge. SPP interference 
causes one contact to ’’pull” light close to its edge and 
the other to ’’push” it away from it. The situation re¬ 
verses at^TOOnm. To confirm that these light ’’attrac¬ 
tion” and ’’expulsion” effects are not related to interfer¬ 
ences within the contact, i.e. that they are independent 
of contact width, we perform the same calculations for 
9, 11 and 13 ridges. While we observe some small inter¬ 
ference effects at longer wavelengths (>730nm, especially 
for d=0), within the primary wavelength range of inter¬ 
est (630-700nm) the absorption is independent of grating 
size. Furthermore, the response is well saturated above 
10 grating periods, consistent with numerical studies on 
the influence of the number of grating periods (4914^. 

The mechanism responsible for these sharp contrasts 
in absorption is interference: the excited SPPs travel 
down the termination step towards the exposed SLG and 
interfere with the incident waves there, as depicted in 
FigE^. Interference between SPP waves was reported 
in pure metal systems (H^. 58, 65, [^. To confirm that 


what shown in FiglTjD is due to interference of SPPs with 
the incident wave, we examine a series of d values, and 
plot the absorption in the exposed SLG within the first 
lOOnm from the contact edge. These calculations are 
performed at 680nm, i.e. in-between the two peaks in 
FiglT]3. An oscillatory response is observed in FigEjD 
(points), which is characteristic of interference between 
two waves with a variable phase between them. The to¬ 
tal field at the SLG will be Einc + Espp, where Fine = 1 
for an incident wave normalized to unit amplitude, and 
Espp = is the SPP amplitude, with A 

the relative SPP electric field strength compared to the 
input field, kspp the SPP wave vector calculated from 
EqlU for the Au/air interface, and a constant phase. 
We thus fit the SLG absorption of FigEjD to: 


A = Aslg 1 + 


(4) 


where jhc = 2.3% is the absorption coefB- where Aslg is evaluated from EqI31 As discussed later, 

dent for suspended SLG in airj^ and ns,02=l-46 is the SPPs on the Au/air interface can leak into the dielectric 
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FIG. 2: (a) Depiction of the interference mechanism, (b) Rel¬ 
ative (i.e. normalized to the case without the grating) SLG 
absorption within lOOnm from the contact edge as a function 
of d. The line is a fit of EqU] which assumes SPP+incident 
wave interference. The top inset color-codes the absorption 
distribution in the exposed SLG as a function of d. (c) Asym¬ 
metric contact layout with d=0 and 410nm, as shown in the 
insets, with GOOnm exposed SLG width. The asymmetric ab¬ 
sorption (i.e. Ar—Al, where Ar, Al are the SLG absorptions 
in the right and left halves of the SLG channel), is propor¬ 
tional to the net photovoltage under coherent uniform illumi¬ 
nation. The insets color-code the absorption profiles for the 
two peaks indicated by the arrows, (d) Asymmetric absorp¬ 
tion when the 410nm contact is replaced by a 770nm one. 


substrate]^ providing an extra loss mechanism. To ac¬ 
count for these losses we scale the imaginary part of the 
wave vector kspp = Re {kspp} ilm {kspp} according 
to Im {kspp} s • Im {kspp}- We treat A, 0o and s as 
adjustable parameters and fit EqH] to the simulation of 
Fig 12^ (note that the oscillation period is dependent only 
on Re {kspp}). An excellent fit (line) is obtained with 
A=1.87, 00 = 7r/5 and 5=9, confirming SPP excitation, 
propagation and interference. FiglJjD (top) color codes 
the absorption within the first 500nm of exposed SLG as 
a function of d. The oscillation of light ” attraction” and 
’’expulsion” from the contact edge is apparent. 

The implication of Eql4] is is that the asymmetry in 
photovoltage is larger than what seen in FiglTjD. E.g., 
if the two contacts shown in FigUlD are placed across 
each other with a 600nm gap, then all the light accumu¬ 
lated from both contacts will be funneled close to only 
one of them. Such a scheme offers great flexibility in 
designing asymmetric contacts suitable for uniform illu¬ 
mination, potentially eliminating the need for different 
metallizations We explore this asymmetric contact 
design in FiglJt (see insets). The asymmetric absorption 
is defined as Aasym = where Ap is the SLG 

absorption in the 300nm area close to the right (d=0) 
contact and Ap is the SLG absorption in the 300nm area 
close to the left (d=410 nm) contact. Aasym is then pro¬ 
portional to the expected net photovoltage. An anti¬ 
symmetric response is obtained, as shown in FiglJt. In 
particular, at 650nm (point (i) in FiglS^) all the flux is 
’’pulled” to the left contact, while at 700nm (point (ii) 
in FigISt) all the flux is ’’pushed” to the right contact. 
I.e., a spectrally-selective region less than 50nm wide is 
created, within which the photovoltage abruptly changes 
sign. Outside this region (A <580nm and A >800nm) 
both contacts have similar responses, thus the net photo¬ 
voltage is zero. Yet, this is not the only response function 
available. Plotted in Figl2]i is the case of a left contact 
with d=770nm. The response is found to be symmetric 
around 660nm. In both cases, the peak absorption is 2- 
3 times higher than what a 600nm SLG on Si02 would 
absorb in the absence of the contacts. In addition, opti¬ 
mizing the contacts just for the highest absorption (i.e. 
without creating ’’clear” symmetric or antisymmetric re¬ 
sponse functions), the peak absorption exceeds 6%, i.e. a 
4-fold increase compared to SLG on Si02. If we further 
limit ourselves to the first lOOnm from the contact edge 
(i.e. within the junction), then Figl2^ gives an 8-fold en¬ 
hancement. Further exploration of the system’s response, 
including absorption in the unexposed (buried) SLG un¬ 
der the contact, different light polarizations, unpatterned 
contacts and Si02 (300nm)/Si substrate, for both uni¬ 
form and focused illumination, are consistent with what 
discussed above (see Methods). 

A detailed analysis of all the SPPs that can be 
launched in our system is performed by considering a 
patterned contact with a large (i=20/im, illuminated by 
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FIG. 3: (a) Field intensity distribution for 650nm focused illumination (l/xm spot size) of a contact with a large (i=20/im 
on one side. (b,c) Field amplitude profiles along the Au-air and Au-Si02 interfaces. (d,e) Spatial Fourier amplitudes of the 
fields along the interfaces, (f) Dispersion relations between frequency and real part of the wavevector. The lines are from Eq[T] 
and the squares are the peaks of the simulated Fourier amplitudes. The vertical and horizontal dashed lines denote the SPP 
condition in Eq[T]with a=620nm. The light lines of the two dielectric media (air and Si02) extend between the vertical axis at 
kx = 0 (normal incidence) to the kx = ncu/c line (90^ incidence), defining the phase-space of allowed free-propagation modes 
in the index n medium 


a narrow l//m width TM-polarized source positioned on 
top of the grating, as shown by the arrow in Figl3^. In 
contrast to plane wave incidence, for which k\\ =0 and 
where only SPPs compatible with kspp = mK are ex¬ 
cited (see Eql2]), the focused beam has |/c||| > 0 incident 
wave vector components. It thus allows the full spec¬ 
trum of SPPs to be excited according to Eql2l Absorb¬ 
ing boundaries are employed in all directions to avoid 
scattered light from re-entering the computational cell. 
In FiglSti the electric-field intensity distribution is plot¬ 
ted at 650nm. Strong scattering and fields extending 
many microns away from the grating are observed, both 
in the Au/air and Au/Si02 interfaces. An intensity oscil¬ 
lation is also observed in the latter, with a period^ 1.4/im. 
SPPs have both longitudinal and transverse field com- 
ponents, the latter being perpendicular to the metal 
surface 45. 4^. Figl3)3,c plots the longitudinal \Ex\ and 
transverse \Ez\ electric field amplitudes at the Au/air 
and Au/Si 02 interfaces, as a function of distance x from 
the grating. A simple decay curve is obtained at the top 
interface, but an oscillating decay curve is obtained for 
the bottom interface. Figl3jl,e report the spatial Fourier 
transform on these fields. A single peak is found in the 
top interface, and two in the bottom one: one at ex¬ 
actly the same wavevector as in the top interface, and 
the other at a larger wavevector. By repeating this proce¬ 
dure at different illumination frequencies we get the SPP 
dispersion shown in FigISf. Lines denote the theoretical 
dispersion curves from EqlU using An’s dielectric func¬ 


tion and assuming either an Au/air or Au/Si02 inter¬ 
face. Squares indicate the peaks obtained from the sim¬ 
ulations by the spatial Fourier transforms. In the Au/air 
interface only the Au/air SPP dispersion curve emerges, 
while both Au/air and Au/Si02 SPP dispersions emerge 
in the Au/Si02 interface. SPPs are surface waves bound 
on a metal/dielectric interface because they exist below 
the light-cone of the dielectric (i.e. the phase-space of 
free propagating modes defined hy k\\ < ncj/c)|45l. l46|. 
The SPP at the Au/Si02 interface is below both the air 
and Si 02 light-cones (i.e. kspp > 

thus cannot couple to any free radiation states. On the 
other hand, the SPPs at the Au/air interface are be¬ 
low the air light-cone and within the Si 02 light-cone 
(i.e. Tiair^jc < kspp < nsi 02 ^!fhus they can leak 
(tunnel) into the free radiation states of substrate, ex¬ 
plaining why we obtain two SPP signals in the Au/Si02 
interface. Such leaky waves have been used for SPP 
characterization within the context of Leakage Radiation 
Microscopy!^ (a far-held optical method analyzing the 
leaked SPP waves in glass substrates to characterize SPP 
propagation on the top interface of a hat or nanostruc- 
tured metal hlm[^). In a semi-inhnite substrate they 
will just propagate away. In a hnite one, on the other 
hand, part of the leaked waves will be rehected back to 
the interface and contribute more to the SLG absorp¬ 
tion. The dominant ehect in photoresponse, however, 
remains in the Au/air SPP, as inferred by the Fourier 
amplitudes in FigslSb-e, and conhrmed by the absorp- 
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FIG. 4: a) Schematic biosensing device. d=600nm is assumed 
for the exposed SLG. An analyte (purple) is deposited on both 
sensor arms, (b) Asymmetric absorption (proportional to the 
net photovoltage) as a function of analyte thickness for four 
sensor arm pairs. The operating wavelength is 678nm. (c) 
Response curve slope as a function of small sensor arm length. 


tion in FiglTjD, with minimal contribution from Au/Si 02 
SPPs at 930nm. 

Besides photodetection, SPP+incident wave interfer¬ 
ence in a MGM architecture also lends itself to label- 
free surface plasmon biosensing ( gsI - ItH . whereby SLG as¬ 
sumes the role of an integrated transducer providing di¬ 
rect electrical readout, thus eliminating the need for op¬ 
tical measurements. The use of SLG as an integrated 
transducer was reported in a dielectric waveguide sensor 
geometry 0, but not in surface plasmon sensing. FigH] 
assumes that one termination step (sensor arm hereafter) 
has length di and the other (i 2 , so that they are at the 
highest slopes of Figl2^, i.e. in the midplane of the 
interference oscillation, with one of them at a positive 
slope and the other at a negative one. In this setup, 
the two contacts are at an accidental degeneracy, pro¬ 




FIG. 5: MGM-PD with plasmonic grating coupler. a) 
Schematic, b) SEM micrograph. 


ducing the same interference between SPP and incident 
wave, thus zero net photovoltage under uniform illumina¬ 
tion. If now the dielectric environment around the sensor 
arms changes by the presence of an analyte, it will cause 
an increase in kspp, thus an additional phase to both 
SPPs. Having the two arms on a different slope in the 
response curve of Figl2^ introduces an asymmetry, thus 
a net photovoltage. The larger the dielectric change, the 
larger the photovoltage. Also, the longer the sensor arms, 
the higher the sensitivity, as the SPP will travel a longer 
distance, therefore sampling more analyte. 

Fig 0] numerically tests this idea simulating four pairs 
of sensor arms, tuned to operate at 678nm according 
to Figdt: (i) (ii=410nm, d2=750nm, (ii) di=775nm, 
d2=1075nm, (hi) (ii=1070nm, (i2=1400nm and (iv) 
di=1435nm, (i2=1730nm. In cases (i)-(iv) the left arm 
alternates from being on a positive to a negative slope of 
the interference diagram of Fig (It. Thus, the photovolt¬ 
age in the presence of the analyte is also expected to alter¬ 
nate. For simplicity, the analyte is assumed to be a thin 
film deposited on the sensor arms and to have an n=1.55, 
an average value for dry protein films[^. FigHjD plots 
the asymmetric absorption as a function of analyte thick¬ 
ness. We obtain a linear response, with slope increasing 
the longer the sensor arms. In FigHt the absolute value 
of the slope is plotted as a function of the small sensor 
arm length di, and a good linear fit is obtained. Tuning 
the arm dimensions thus provides an additional tool for 
controlling and tuning the device’s performance and sen¬ 
sitivity. At long arm lengths, SPP losses will limit the 
sensitivity, but they could be overcome, e.g. by increas¬ 
ing the metal thickness and/or reducing the substrate 
index to limit SPP coupling to the substrate. 

Having demonstrated the design versatility of grating- 
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FIG. 6: Scanning photovoltage maps for excitations at a) 514nm, b) 633nm, c) 785nm and perpendicular polarization (TM). 
The grating coupler dimensions are not to scale in the vertical direction, d) Photovoltage at 785nm as a function of polarization. 
0° denotes TM polarization. 


coupled GPDs, we now turn to the experimental valida¬ 
tion of our predicted SPP enhanced SLG absorption and 
photodetection. The general architecture of our experi¬ 
mental devices with plasmonic grating coupler is shown 
in Figl5^. SLG is contacted with metallic source and 
drain electrodes. In order to quantify the enhancement 
relative to the no-grating case, we fabricate the grating 
coupler on top of one of the contacts and leave the other 
contact flat. This is the simplest way to break the con¬ 
tacts symmetry and allows the generation of net non-zero 
photoresponse, even with both contacts illuminated. 


Our devices are fabricated as follows. Graphene is 
produced by mechanical exfoliation of graphite onto 
an Si-T300nm Si02 3, ^ and characterized by optical 
microscopy [Tij and Raman spectroscopy [ t^ Sub¬ 

sequently, the source and drain contacts are prepared 
by e-beam lithography and a base metallization layer is 
deposited by thermal evaporation of 4nm Cr and 50nm 
Au, employing a lift-off step. The 620nm period grating 
is then defined in a further e-beam lithography step by 


performing a second thermal evaporation of 50nm Au fol¬ 
lowed by lift-off. FigISb shows a scanning electron micro¬ 
graph (SEM) of the device. A slight asymmetry between 
the ridges and grooves is detected due to overexposure 
during e-beam lithography, but this does not change the 
spectral characteristics of the grating, solely determined 
by its period. Afterwards, the samples are bonded into a 
chip carrier for electrical and optical characterization. 

We perform wavelength dependent photovoltage map¬ 
ping to determine the spatial pattern of the devices’ pho¬ 
toresponse. FigE] plots the photovoltage maps at differ¬ 
ent incident wavelengths for polarization perpendicular 
to the grating (TM-polarization, lOOx ultralong working 
distance objective, numerical aperture NA=0.6). FigE] 
also shows the structured grating contact and the fiat 
contact without perturbations. At 514nm the photore¬ 
sponse occurs predominantly at the contacts’ edges and 
is of similar magnitude, but opposite polarity, to that 
of a standard MGM PD[^, Q. At 633nm, FigEb in¬ 
dicates that the influence of the grating starts emerg- 
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(a) 


illumination 



exposed 


SLG 

4nm 

Si02 


Si 



FIG. 7: Calculated absorption in the SLG exposed part when 
a l/xm spot size illuminates the grating 4/im away from the 
metal/SLG junction, (a) Schematic of the simulation system. 
300nm Si02 on Si is assumed for closer matching to the ex¬ 
periments. (b) Absorption in SLG. Two polarizations and a 
flat contact case are considered. 


ing. 2-3/im away from the edge of the patterned contact, 
a photoresponse is visible, even though no junction is 
present. The effect is much more pronounced at 785nm, 
where the entire structured contact becomes photosensi¬ 
tive, Figinb, and the photoresponse is enhanced^400% 
compared to the flat contact edge. Furthermore, the 
responsivity of the device is polarization dependent, as 
shown in Figibji. The strongest photoresponse occurs for 
perpendicular polarization (0°,TM-polarized light). 

To have a direct comparison with the experiments, 
we perform calculations with a focused beam illumina¬ 
tion (lyum width) on top of the grating,^4//m away from 
the contact edge, as depicted in FigC^. The termina¬ 
tion step length is taken as 1250nm (estimated from the 
SEM image in FigEJ and we also include the Si sub¬ 
strate with 300nm Si 02 . FigI 33 plots the absorption in 
the exposed SLG 4//m away from the illumination spot 
for TM polarization, TE polarization and a flat contact. 
Both the theoretical spectral and polarization responses 
are in excellent agreement with experiments, and ver¬ 
ify the strong responsivity above 700nm. This contrasts 
with the normally incident plane wave illumination case, 
below 700nm (see Eig(T] 3 ). This is understood by con¬ 
sidering Eql2]and Eigl3f. Eor a normally incident plane 


wave and focused illumination, k\\ is zero and nonzero, 
respectively. Since the latter case is less restrictive for 
SPP excitation, it results into a wider SPP spectrum at 
both Au/air and Au/Si 02 interfaces, thus into a wider 
responsivity compared to the plane wave case (see Meth¬ 
ods for details) 

In conclusion, we demonstrated the coupling to 
graphene of surface plasmon polaritons excited in 
a metallic plasmonic grating and its exploitation in 
graphene-based photodetection with enhanced respon¬ 
sivity and polarization selectivity. Depending on its 
dimensions, highly tunable spectral selectivity below 
50nm bandwidth can be achieved. Further, the sym¬ 
metry of the photodetector can be broken making it 
operable under full illumination, despite identical metal 
source and drain contacts. The underlying mechanism 
involves the coupling of light into SPPs on the patterned 
contact, and their propagation to the exposed SLG 
area. For uniform coherent illumination, these SPPs can 
further interfere with the waves directly incident on the 
exposed SLG, offering a novel tuning capability where 
the light flux can be attracted or repelled from the 
contact edge by design. The whole contact thus becomes 
a highly tunable polarization- and spectral-selective 
photosensitive area. SPPs and incident wave interfer¬ 
ence can potentially be employed for (bio-) sensing by 
tailoring the grating dimensions. This may allow a novel 
plasmonic sensing architecture with high sensitivity 
and small footprint with direct electrical readout and 
without complicated optics. 

We acknowledge funding from EU Graphene Flagship 
(no. 604391), ERG Grant Hetero 2 d, EPSRC Grants 
EP/K01711X/1, EP/K017144/1, EU grant GENIUS, a 
Royal Society Wolfson Research Merit Award. 


METHODS 


FDTD simulations 


In our FDTD simulations Maxwell’s equations are 
time-integrated on a computational grid, with a Drude- 
Lorentz model assumed for the dielectric function]^: 


e{uj) = € 00 - 


UJ" 


N 




— 7 


• 1 

j = l 3 


iujV ^ 


(5) 


where the first term is the Drude free-electron con¬ 
tribution and the second contains Lorentz oscillators 
corresponding to interband transitions, ujp and I /7 
are the free electron plasma frequency and relaxation 
time, Uj, Aej, F^ are transition frequency, oscillator 
strength and decay rate for the Lorentz terms. To 
accurately reproduce the experimental dielectric func¬ 
tions (An from Ref.j^ and Si from Ref. 0 ) we treat 
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FIG. 8: Refractive indexes used in the calculations for a) An, 
b) Si. Symbols are experimental data from Refs.[^.[^: lines 
are the corresponding Drude-Lorentz fits. 


these as fit parameters. For Au we use A/'=4, and 
e^=4.054, Ae^=(0.43, 0.634, 0.755, 1.059), hujp=8.76eV, 
hj=0m8eV, hVLj={2m, 3 . 03 , 3 . 54 , 4.23)eV, and 
nF^.=(0.458, 0.641, 0.892, 0.959)eV. For Si we use N=1 
without any Drude term: eoo=l-89, Aej = (1.198, 0.963, 
1.021, 1.164, 1.407, 2.259, 1.869), fifl^=(3.39, 3.51, 3.68, 
3.86, 4.06, 4.25, 4.61)eV, and ;iFj=(0.188, 0.203, 0.239, 
0.269, 0.283, 0.265, 0.0)eV. Figl8] plots our model dielec¬ 
tric functions along with the experimental ones, showing 
an excellent agreement. 

The computational cell comprises a 2nm cubic grid. 
Perfectly matched layer (PML) boundary conditions (soj 
are applied at the edges of the cell in the vertical direc¬ 
tion. In the lateral direction, we use periodic boundary 
conditions for the uniform illumination case, and PML 
for the focused. 


SLG index of refraction 

The index of refraction of SLG was measured by ellip- 
sometry for wavelengths up to 750nm in Ref.j^. To ex¬ 
trapolate to longer wavelengths, we assume undoped SLG 
and use the universal optical conductance a = e^/4fi, to 
get the dielectric function: 

. 47rcr . aX 

e = eoo+^—-= Coo+^773- (6) 

^CtSLG ^CiSLG 

where (i5'LG=0.335nm is SLG’s thickness and a = e‘^/hc 
is the fine structure constant. In the thin film limit, Eql6] 
yields: 


= Im(e) = na = 2.3% (7) 

A value of eoo = 5.7 ensures that EqE] matches the 
ellipsometric experimental data of Ref.[63| at smaller 
wavelengths. We fit a Drude-Lorentz model to both 
the ellipsometric data at small wavelengths and EqE] 




FIG. 9: (a) SLG refractive index used in our calculations. 
Symbols are ellipsometric data from Ref. . (b) SLG ab¬ 
sorption corresponding to the applied model. 


at longer wavelengths, Figl9^. In Fig|9]3 the resulting 
SLG absorption is plotted. The Drude-Lorentz model for 
SLG uses N=4 with eoo=2.148, Ae^=(64.8, 2.92, 1.69), 
fico’p=1.34eV, fi7=0.7eV, hQj={1.0^ 4.0, 4.56)eV, and 
fiFj = (5.41, 2.77, 1.0)eV. Within FDTD’s 2nm grid, SLG 
is treated as an effective 2nm thick slab, thus its dielectric 
function is scaled to reproduce the correct absorption and 
reflection properties according to e ^ 1 ^ {e — l)dsLG/‘^^ 
i.e. €qq —^1.192 and Ae^ —^(10.85, 0.488, 0.283). 


SLG absorption in the thin film limit 


In a three-layer system, the normal incidence Fres¬ 
nel ^nations for reflection and transmission amplitudes 
are 81, 82| : 


r 


t 


1 - ri2r23e^*'^ 


^12^236^"^ 

1 - ri2r23e^'^ 


( 8 ) 

(9) 


where Vij = — rij)!{rii + n^), Uj = 2ni/{rii -h nj ) and 
(f) = n 2 Ujd 2 lc^ with ^2 the film thickness. Here we assume 
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FIG. 10: Absorption in the l/xm exposed part of SLG (red line) and in the covered SLG under the An contact (black line). 
The green dotted line is the absorption in l/xm of SLG on top of Si02. All curves are normalized to the light flux incident 
in l/xm. (a) TM-polarized, (b) TE-polarized. (c,d) The two polarizations for unpatterned contacts. In cases (b,d) a peak 
at~500nm is observed. This is related to a transparency window in An, as shown in the inset where the reflection, absorption 
and transmission through a 50nm thick An film suspended in air are plotted. 


the incoming medium ni = Uair = 1, the film 77-2 = nsLG 
and the semi-infinite substrate ns = nsi 02 ' Reflection 
and transmission coefficients are R= |rp and T = nsjtp. 
In the thin film limit, using the SLG dielectric function 
given by EqEl we get: 


\ni + n 3 / \ nf-n§ 


/3(e - nina) ^ V 
9e(n3 - ni) y 


4nin3 f _ 2/3 \ 

(ni + nzY V ni+nsj 


( 11 ) 


Regarding reflection, this is typically dominated by the 
substrate since the difference in refractive index is usually 
largest between air and the substrate. However, in the 
ni = ns = n case the third term in the right hand side of 
EqUOjis the only nonzero term, resulting into NLG reflec¬ 
tion: R = — n^p/(n^A^). This is very small 

for SLG {R = 0.02% at 600nm for suspended SLG [13) 
but it increases quadratically with N. 


Coherent Uniform illumination 


where /3 = 27rdlm{e} /A = Nna, e is the film dielectric 
constant, a = e^/he is the fine structure constant, and 
d = NdsLG is the N-layer graphene (NLG) thickness. 
This results into the absorption: 

A=l-R-T^NTia- -^ (12) 

{ni^nsY 


while for a NLG suspended inside a uniform medium (i.e. 
ni = ns = n) we find A = Nnajn = 2.3/n%. These 


equations are valid for Nil0[27|, l76|, l83|. For > 10 the 
thin film limit breaks down and the optical paths inside 
the film must be taken into account. 


Fig do] plots the response of our system considering the 
absorption in both exposed and covered SLG for both po¬ 
larizations and patterned and unpatterned contacts. For 
simplicity we only consider the (i=410nm case with 13 
ridge periods. We find that in all cases, other than TM- 
excitation on the patterned contact, the absorption in the 
exposed SLG is similar to that of a SLG on top of Si02 
in the absence of the contact (green dotted line). Inter¬ 
estingly, there is some absorption in the covered SLG as 
well. In the TM-excitation of the patterned contact we 
obtain a similar modulation with wavelength as that in 
the exposed SLG, albeit of significantly smaller magni- 
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FIG. 11: (a) Field intensity distribution for a l/xm spot il¬ 
lumination at normal incidence (TM), for 500, 600, 700nm. 
The color-coding is in logarithmic scale, (b) SLG absorption 
distribution (vertical axis) as a function of illumination wave¬ 
length (horizontal axis). The three wavelengths studied in 
(a) correspond to the three characteristic cases: absorption 
under the illumination spot, extended absorption under the 
contact, and absorption in the exposed SLG far from the illu¬ 
mination spot, (c) the same as in (b) but for TE polarization, 
(d) Same as above for unpatterned contacts (both polariza¬ 
tions yield identical results). Absorptions are normalized to 
the incident light flux. 


0.6 
0.4 
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0.0 
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FIG. 12: (a) Absorption in the exposed SLG for the system 
studied in Fig.l with a SiO2(300nm)/Si substrate, (b) Ab¬ 
sorption interference enhancement without An contacts. The 
modulation at~570nm is due to I-E. 

tude. We note that, in order to directly compare with the 
exposed SLG absorption, we still normalize to the power 
illuminating the exposed SLG area. Had we normalized 
to the flux illuminating the whole contact, then the cov¬ 
ered SLG absorption would appear with a much smaller 
magnitude than in FigJTOl We also obtain absorption in 
the other cases as well, peaked at^SOOnm. This, however, 
is irrelevant to the grating and SPPs, as it is there that 
An becomes most transparent. An is strongly absorptive 
at high energies because of the onset of interband transi¬ 
tions from its d-electrons[^ while it is strongly reflective 
at small energies because of its conduction electrons}^. 
The inset in FigJTOtl plots the reflection, transmission and 
absorption coefficients through a 50nm thick An film in 
air, displaying An’s transparency window at^bOOnm. 

Focused illumination 


Wavelength (nm) 



J_1_I_I_L 


We explore the case of focused illumination on the grat¬ 
ing, which more closely resembles the experiments. To 
better facilitate the simulation and avoid having laterally 
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FIG. 13: System of Fig.3 in with a Si02 (300 nm)/Si substrate. Besides a small modulation in the Fourier amplitudes, the 
SPP characteristics remain the same as for the semi-infinite Si02 substrate. Notable is also a deeper ’’beat” modulation of the 
field intensity at the Au-Si02 interface, because of reflection of the ’’leaked” SPP at the Si02/Si interface. 


scattered light re-enter the computational cell, we remove 
the lateral periodic boundary conditions and adopt PML 
boundary conditions [s^, so that any light scattered to¬ 
wards the sides of the computational cell permanently 
exits the calculations. We also consider a much larger 
exposed SLG area for better visualization. We adopt a 
TM polarized l/rm-wide plane source illuminating the 13- 
period (i=410nm grating at a non-symmetric position, as 
depicted in FiglTTl The frequency domain electric-field 
intensity profiles for A=500, 600 and 700nm are shown 
in logarithmic scale in FiglTTk. No scattering occurs at 
500nm, while the most intense is seen at 700nm. 

The system’s full response is shown in FigITTb. where 
we plot the SLG absorption throughout the length of 
the structure (vertical axis) for different wavelengths 
(horizontal axis). The illumination source is again a 
TM-polarized Ijnm wide spot as shown in the inset 
schematic. Absorption is normalized to the peak incom¬ 
ing flux per unit area. Three distinct regions emerge: 
up to 550nm, absorption only occurs in the covered 
SLG directly underneath the illumination source. For 
550nm< A <650nm, increased absorption is found in an 
extended area (several microns away from the illumina¬ 
tion spot) in the covered SLG, pointing towards light 
diffraction into SPPs in the Au-Si02 interface. Above 
650nm there is strong absorption in the exposed SLG 
beyond the grating. This points towards light diffraction 
into SPPs in the Au-air interface, which propagates and 
reaches the exposed SLG at the grating’s edge. In FigITTb 
we plot the same absorption map for TE illumination of 
a patterned contact, while in FigITTb for illumination of 
an un-patterned contact. In both these cases no absorp¬ 


tion is found, except directly underneath the illumination 
spot. It is thus clear that SPP-mediated effects are dom¬ 
inating the response for TM-polarized illumination. 


Si 02 /Si substrate effects 

The Si02 (300nm)/Si substrate has two effects. First, 
it provides some interference-based enhancement in the 
SLG absorption. FiglT^ plots the exposed SLG 
absorption for the system described in Fig.l on top of 
a 300nm Si02/Si substrate. The response is similar to 
the semi-infinite Si 02 case, except for an overall modula¬ 
tion due to the interference effects in the Si 02 dielectric 
spacer. The net interference-enhancement (I-E) effect on 
absorption (i.e. without the patterned contact) is plot¬ 
ted in EigHSb). The dielectric spacer can thus be used as 
an additional degree of freedom (i.e. the spacer’s index 
and thickness) in optimizing the system’s response. We 
also simulated the asymmetric contact layouts studied 
in Eig.ld on top of 300nm Si02/Si and found that they 
produce a very similar response irrespective of the finite 
dielectric spacer. This is expected, since I-E partially 
cancels out when considering asymmetric absorption. 

The second effect of the 300nm Si02/Si substrate is 
that the leaked SPP of the Au-air interface will be re¬ 
flected back from the Si substrate. EigUd] plots the field 
intensities along the two An interfaces as well as the cor¬ 
responding Eourier transform amplitudes (see Eig.3b-e 
for the 300nm Si02/Si case). We note a deeper ’’beat” 
modulation of the fields on the Au-Si 02 interface, due 
to the back reflected fields of the ’’leaked” Au-air SPP. 
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This is also apparent in the Fourier amplitude. However, 
no frequency shifts are observed in the latter. We thus 
conclude that, other than some small amplitude modula¬ 
tions, the SPP structure is largely unaffected by the Si02 
being 300nm thick or semi-infinite. 


* Electronic address: acf26@eng.cam.ac.uk 
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